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ABSTRACT 

An algorithm has been developed that converts penetration angles reported by Static Earth 
Sensors (SESs) into Earth observation vectors. This algorithm allows compensation for variation 
in the horizon height including that caused by Earth oblateness. It also allows pitch and roll to 
be computed using any number (>1) of simultaneous sensor penetration angles simplifying 
processing during periods of Sim and Moon interference. 

The algorithm computes body frame unit vectors through each SES “cluster”. It also computes 
GCI vectors from the spacecraft to the position on the Earth’s limb where each “cluster” detects 
the Earth’s limb. These body frame vectors are used as sensor observation vectors and the GCI 
vectors are used as reference vectors in an attitude solution. The attitude, with the unobservable 
yaw discarded, is iteratively refined to provide the Earth observation vector solution. 

INTRODUCTION 

Static Earth Sensors 1 (SES) can be used to determine the pitch and roll attitude of Earth pointing 
spacecraft by detection of the amount of light entering several photocells (eyes). These eyes are 
directed so that at nominal altitude and attitude the Earth edge intersects with their fields-of-view 
(FOVs). Any attitude motion will move the eyes with respect to the horizon so that the portion 
of the FOV intersecting the Earth will change. Therefore the total light (and therefore output 
current) of the various eyes is a function of the attitude. Because spacecraft yaw moves each eye 
around the horizon, the SES is quite insensitive to yaw. 

A simplified description of the SES serves to illustrate the operation of an SES and provides a 
definition for some of the terms that will be used later. SESs generally contain clusters of 
photocells. For complete observability, at least two clusters are required, and the vectors from 
the body origin to these two clusters can not be coplanar with the Z-axis. 

Each cluster contains one or more eyes, usually with triangular FOVs. These eyes measure the 
position of the Earth limb relative to the spacecraft and, when there are more than one eye, such 
quantities as the background light level. SESs typically have a passband in the infrared (14-16 
microns). In this region of the spectrum the Earth brightness is almost independent of the 
presence or absence of solar illumination. 

The configuration of ak^ES containing four clusters with four eyes each is shown in Figure 1. 
The four eyes in each cluster include one that views space and provides a background level, one 
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Figure 2. View From Spacecraft to Earth for a Typical SES 


that is the primary Earth detection eye (with its broad base towards the Earth) and two that 
enhance the response of the sensor. The eyes have sides that are typically 5 deg in extent. 

The amount of light entering the eyes is expressed as a penetration angle, 8 . The penetration 
angle can be thought of as the angle between the edge of the Earth image in the FOV of the 
central eye and a reference point. This is illustrated for a simplified SES in Figure 2. 



Figure 2. Definition of Penetration Angle For SES 
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For each cluster, C, define a sensor coordinate system: X c , Y c , and Z c so that: 

Xc=~ 


\Jt c xZ\ 

and 

Z C =Z 


ZxR c 

ZxR c 

X c xZ 
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where Z is the body Z-axis and R c is a vector from the body origin through the reference point 
of the cluster. The sensitive direction of the sensor is in the Yc-Z plane. 

The penetration angle for cluster C is almost exactly proportional to rotation about X c (the 
sensitive axis for this cluster) and is relatively insensitive to rotation about Y c (the non-sensitive 

axis for this cluster). The approximations that the penetration angle is proportional to the first 
rotation and independent of the second limit the accuracy of the algorithm described in this 
paper. 

For an SES with four clusters, spaced uniformly around the spacecraft Z-axis (the axis that is 
nominally directed towards Earth), clusters on opposite sides of the spacecraft have antiparallel 
Y c axes while adjacent clusters have X 2 = Y x (with equivalent relationships for the other 

clusters). Defining rotations about 7, as 6, and those about Y 2 as <f>, the attitude can be computed 
from the penetration angles by: 


2 

( 4 ) 

2 

( 5 ) 


where the subscripts indicate the cluster measuring the penetration angle. 6 and <j> can be easily 
converted to body frame roll and pitch angles using the relationship of the SES coordinate 
system to the body. Equations (4) and (5) and are used under normal conditions for spacecraft 
such as the Tropical Rainfall Measurement Mission 2,3 (TRMM). 

There are two difficulties with this simple approach to SES attitude determination: compensation 
for horizon height variation and generality. 

Eqs. (4) and (5) are based on the comparison of angles from the Z axis to the detected horizon 
with angles from the nadir to the horizon. If the Earth appeared spherical in the infrared, nadir- 
horizon angles would be independent of the location on the Earth where the horizon was 
detected. The Earth’s oblateness causes a latitude and yaw dependent variation in these angles. 
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Seasonal and position dependent variations in stratospheric temperatures also cause variations of 
nadir-horizon angles but, unlike oblateness effects, are only approximately predictable. 4 

Oblateness and (if an accurate model is available) stratospheric temperature effects can be 
compensated by proper adjustment of the penetration angles based on the position and 
orientation of the spacecraft. This compensation, however, is complex and mission dependent. 

Eqs. (4) and (5) are applicable to a specific sensor configuration. The SES must include four 
clusters and the clusters must be mounted at 90 deg. intervals around the Z-axis. In principal, 
only two clusters are needed to obtain roll and pitch information and any greater number could 
be used. In practice, even with four clusters mounted properly, there are times when one or more 
of the clusters is unavailable (due to Sun or Moon interference). During these times, Eqs (4) and 
(5) must be replaced by more complex algori thms 5 which are generally also specific to this 
sensor geometry. The Aqua spacecraft, for example, has an SES with two, non-orthogonal 
clusters, and consequently requires a different algorithm from these equations. 6 

This paper describes an algorithm that can be used to determine attitudes for a general SES. It is 
independent of the number and placement of clusters. Using the algorithm, compensation for 
deviations of the Earth’s shape from a sphere is relatively simple. It requires only a model 
providing the effective Earth radius as a function of time and location on the Earth’s surface. 

ALGORITHM 


The algorithm to determine pitch and roll from SES measurements is iterative. Each iteration 
consists of the following procedure: 


1. At each time, convert measurement data for each SES cluster into body frame horizon 
observation vectors, O c . These vectors go from the body center through the point in each 
eye at which the Earth horizon is detected. 

The body frame vectors, R c , are known from the sensor mounting geometry. The body 


frame observation vectors are calculated by rotation of these vectors around the X c -axis by 
the penetration angle: 

OcjBodv - R°t(Xc.Body > $c )^C ( 6 ) 

2. At each time compute Geocentric Inertial (GCI) reference vectors, H c ,from the spacecraft 

body to the point on the horizon at which the Earth would be expected to be detected using 
current estimates of the attitude. 

2. 1 . Initially, use nominal values for the pitch and roll. For yaw, use either nominal values or 
values obtained from another sensor. On subsequent iterations use the most recent 
attitude estimates for pitch and roll. 

2.2. Using the attitude, convert the vectors Xc,Body to GCI: 


X-cfici ~ A X CBody 


( 7 ) 


2.3. Compute a horizon crossing vector, H caa , nearest to Y c from the intersection of two 
cones. The two cones are: 

one centered on X C Ga with a 90 deg half-angle 


4 


one centered on the nadir vector, with-half-angle, p, corresponding to the Earth 
width as viewed from the spacecraft: 


p = sin 


f r r +Ah'' 


( 8 ) 


where: 

- r E is the Earth equatorial radius, 

- Ah is a deviation from the equatorial radius due to the nominal height in the 
atmosphere at which the sensor detects the Earth, the Earth’s oblateness, and 
any corrections for stratospheric temperature effects 

- r s is the length of the spacecraft nadir vector, N , to the center of the Earth. 

2.4. Scale the horizon crossing unit vector to its appropriate length by: 

Hc.gci = r s c °s(/?)i/ C Ga (9) 

2.5. Compute the vector from the Earth center to the position on the Earth at which the 
horizon crossing occurs by: 

Pc.gci=H c -N ( 10 ) 


2.6. Obtain Ahc from P c and some horizon height model. The horizon height may be 

constant, or variable, based on the time, latitude, and possibly longitude of each horizon 
crossing position. Normally, it will include oblateness compensation at a minimum. A 
separate Ahc is computed for each observation. 

2.7. Compute new values of pc and H c (using the values of Ahc for Ah in Eq. (8)), by 

repeating steps 2.3 and 2.4. Only a single repetition of these steps is needed because Ah 
does not vary greatly with position (see below). 

3. Solve for attitude using all valid H c and O c vectors. The H c vectors are treated as 
reference vectors and theO c vectors as observations. The attitude is computed using the 
QUEST 7 algorithm. 

4. Since the yaw rotation is very poorly observable by the SES alone, rotate the computed 
attitude so that its Z-component is nominal or equal to the input yaw. 

5. Repeat steps 2 to 4 until the attitudes no longer change. This generally requires 3 or 4 
iterations. 

6. Report the X and Y rotations of the final attitude as roll and pitch. 

RESULTS 

The algorithm was tested using simulated data. Two arrangements of clusters were simulated. 
The first arrangement has four clusters at 90 deg intervals around the Z-axis. They are located 
equidistant from the X and Y axes. In cases where only two clusters were used, adjacent clusters 
were selected. For most of the tests, the two clusters were a pair on the same side of the Y-Z 
plane, but in one test the two chosen were on the same side of the X-Z plane. 

The second arrangement for testing used three clusters. This arrangement placed one cluster in 
the Y-Z plane and the other two at 120 deg intervals around the Z-axis. 

The simulation used an ephemeris from the Aqua mission. Aqua’s orbit is nearly circular with a 
semi-major axis of 7070 km and an inclination of 98.2 deg. 
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Simulated truth attitudes were created from ten-term Fourier series using random coefficients. 
Separate coefficients were used for roll, pitch, and yaw attitude components. The range of the 
coefficients was chosen so that the minimum period of any term was 18 deg and the total 
amplitudes of the attitudes were about 1 deg. 

The pitch, roll, and yaw angles created in this way had maximum values of 1.10, 1.02, and 1.07 
deg. and minimum values of -1.24, -1.11, and -1.19 deg, respectively. Their standard deviations 
were 0.404, 0.352, and 0.402 deg. 

The simulated attitudes were sampled at 4 second intervals over a 24 hour period providing 
21601 points for each test case. For each point, the penetration angles were simulated using the 
Aqua ephemeris and an oblate Earth model with a constant horizon height of 30 km above the 
oblate spheroid. An example of the truth attitude for a portion of the data is given in Figure 3. 

In addition to the number and geometry of the clusters, three parameters were varied: 

• The number of iterations of the entire algorithm was varied to determine convergence. The 
number of iterations for finding the horizon height at each location (step 2) was kept at 2. 
Separate tests showed that the maximum change in Ah was 20.82 km for the first iteration of 
step 2, 0.095 km for the second, and 0.006 km for the third. 

• In some cases, a yaw value of zero was used in the simulation while in others the variable 
yaw (described above) was applied. 

• In some cases a zero-mean, random, gaussian noise with standard deviation of 0.02 deg was 
applied to the penetration angles. 

Examples of the results of some of the cases studied are shown in Figures 4-7 and a summary, 
providing the maximum absolute differences between truth and calculated attitude components, 
as well as their standard deviations and the t im e taken for the computation, is given in Table 1 . 
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Figure 4. SES Attitude Errors: Case with 4 Clusters, Nominal Yaw, and No Noise 
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Figure 6. SES Attitude Errors: Case with 4 Clusters, Simulated Unknown Yaw, and 0.02 deg (Ict) Noise 



Figure 7 SES Attitude Errors: Case with 2 Clusters, Simulated Unknown Yaw, and 0.02 deg (Is) Noise 


8 










Table 1. Results of SES Attitude Computations 
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* 2 Clusters on the same side of the X-Z plane (for ail other cases with 2 clusters, they were on 
the same side of the Y-Z plane) 


Several important results can be seen in Table 1. 

1. The algorithm converges after 3 or 4 iterations (In cases 1-6, and 7-12 the errors reach a 
minimum after 3 or 4 iterations). 

2. The inherent l-o accuracy of the method is better than 0.0002 deg for four clusters (case 4), 
about 0.0004 deg for three clusters (case 20), and about 0.006 deg for 2 clusters (case 9). 

3. The geometry of the clusters greatly affects the accuracy — for the geometries used in the 
simulation, and using 2 clusters, calculated values of rotations about an axis between the 
clusters were an order of magnitude more accurate than those about an axis perpendicular to 
the cluster bisector (cases 10 and 13). 

4. If the (unobservable) yaw is not known from other sensors and is not small (« 1 deg) the 
accuracy decreases to about 0.003 deg for the four (case 16) and three cluster (case 22) cases, 
and to as much as 0.007 deg for the 2 cluster case (case 17). 
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5. Addition of noise (0.02 deg la) produces errors equal to the noise for the two cluster case 
(case 15 and 19), and approximately equal to the noise divided by the square root of 2 for the 
four cluster case (case 14 and 1 8). For the three cluster case introduction of noise caused 
errors only slightly greater than those for the four cluster case (compare cases 21 and 23 with 
cases 14 and 18). 

6. The noise is the major overall error source (compare cases 14, 16 and 18 for four clusters, 
cases 15, 17, and 19 for two clusters, and cases 21, 22, and 23 for 3 clusters). 

7. The current implementation of the algorithm can process about 200 observations per second 
on a personal computer operating at 833 MHz. The most time consuming portions of the 
code are retrieving the ephemeris and solving for the attitude with the QUEST algorithm. 

Even with no applied noise and perfect knowledge of yaw, the computed pitch and roll values 
contain small, but noticeable errors. The following discussion shows why the approximation that 
the roll and pitch measurements are independent yield small, but finite, errors. For example, if a 
cluster is in the X-Z plane, its penetration angle is assumed to be independent of rotations about 
the X-axis and proportional to rotations about the Y-axis. The error in this approximation is 
illustrated in Figure 8. 

Figure 8 shows the Earth occultation of the central eye in an SES cluster in the body X-Z plane 
such that Y C =X. This cluster is sensitive to rotations about the Y-axis and should be 
insensitive to rotations about the X-axis. Two cases are shown: in the FOV on the left, the 
X-axis rotation (<J>) is zero, while on the right it is significant. The point A is on the X-axis while 
B and C are on the Earth limb, at the center of the FOV, for the two cases. Note that all of the 
angles have been exaggerated to illustrate the effect. Measured penetration angles are 
proportional to the total occulted area (shaded in the figure). In this algorithm, modeled 
penetration angles are the angles from the reference point (on the base of the triangular eye — see 
Fig 2) to the intersection of the Earth limb with the normal to the base of the FOV (points C and 
B for the two cases). 


When <5 is zero, rotations about Y move the FOV so that the Earth penetrates it to a degree 
proportional to the rotation angle. For ® not equal to zero, two effects occur: 1) the position of 
the Earth limb intersection changes (points B and C) so that for the same Y-rotation a different 
penetration angle is observed; 2) the shape (and area) of the Earth occultation area changes. 

For near no mina l attitudes, the angle from the body X-axis to Nadir (not shown) is nearly 90 deg. 
The angle (<}>) from nadir to the Earth limb at B.is nearly independent of O so the angle (X) from 
the X-axis to B is given by: 


sin/t = 


cos^ 

cos® 


(ID 


If the X-axis is 90 deg from nadir, a rotation about it of 1.5 deg (e.g., <E> = 1.5 deg) results in a 
change of X, and therefore the penetration angle, of about 0.01 deg. 
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Figure 8. Effect of Rotation on Penetration Angle 


The second effect is illustrated in Figure 9. This figure illustrates a trapezoidal portion of an 
FOV with internal angle a. The Earth limb (approximated as straight dashed lines) is drawn for 
a nominal case (line BE) and a case with a rotation of <X> (line AD). The detected penetration 
angle is proportional to the area under BE in the one case and under AD in the other. Using 
plane geometry and the tangent to the Earth limb, instead of the limb itself, as close 
approximations, the area of triangle ABC is: 


AREA(ABC) = 


r 2 sin(<t>)sin(a) 
2sin(a + 0) 


( 12 ) 


where r is the length of CB (= length of CE). 1 ne area of triangle CDE is: 


AREA(CDE) = 


r 2 sin(<l>)sin(a) 
2sin(<ar -<t>) 


The occulted area is: 


AREA(occulted) = 




tan(a) 


(13) 


(14) 


where b is the length of the base of the eye. 

For typical geometry parameters (a = 60 deg, r = 2 deg, b = 5 deg), and deg, the difference 
of the triangles’ areas is on the order of 0.01 percent of the occulted area. This would produce an 
error of the same size in the measured penetration angle. 

Reference 3 contains a more complete description of these approximations. 
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Figure 9. Illustration of Change in Detected Earth Area With Rotation About Non-sensitive Axis 


CONCLUSIONS 

Although the algorithm presented here does not provide exact results, the error inherent in the 
algorithm is smaller than SES measurement errors. The largest component of SES error comes 
from the variation of the Earth limb height. 8 Although part of this variation is systematic, 
varying with latitude and season, the stochastic, correlated portion can change the Earth limh 
height by ±5 km or more. 9,10 At an altitude of 700 km, this results in an unpredictable variation 
in penetration angles of about ±0.1 deg. If such errors were distributed uniformly, their standard 
deviation would be about 0.06 deg. 

In this study the simulated noise led to penetration angle errors of 0.02 deg, about 1/3 of the 
typical values at this altitude. Even at this low level, these variations were the principal source 
of error. In addition, for a well controlled spacecraft, the pitch and roll would be significantly 
closer to nominal, and the yaw known more accurately, than was modeled in this study. These 
factors would further reduce the errors due to the algorithm, as compared to those due to 
unpredictable variations. 

The major advantage of the new algorithm is its generality. 

The same algorithm can be used for Static Earth Sensors, regardless of the number or placement 
of the clusters. No special modification of die algorithm is needed if clusters fail or during 
periods of Sun or Moon interference. 

Several refinements and related applications could be developed. 

Allowing variable cluster weights might improve the method’s accuracy. The weights could be 
chosen based on predicted horizon height uncertainties at the time and location of each Earth 
horizon measurement. 

A variation on this algorithm could be developed for scanning Earth sensors. Scanning Earth 
sensors view the Earth limb by rotating an infrared sensitive detector around an axis so that its 
FOV intersects the Earth at two points. The position of these two points in the body frame can 
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be computed from raw sensor data and treated in the same way that the observations from the 
SES clusters are in the method described here. 

It would also be possible to use horizon crossing vectors directly in an extended Kalman filter. 
This would eliminate the conventional intermediate step of first computing observed nadir 
vectors. In order to use horizon crossing vectors directly, it would be necessary for the filter to 
take account of the different uncertainties in the vector components. 
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